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Background: Repetitive transcranial magnetic stimulation (rTMS) of the dorsolateral prefrontal 
cortex (DLPFC) appears to have short-term antidepressant properties. The aim of the current 
study was to update our previous meta-analysis and to investigate factors associated with the 
antidepressant properties of rTMS. 

Method: Following a systematic literature search conducted in Medline and Psyclnfo, N=14 
sham-controlled, parallel design studies (published after 2008 to August 2013) that had utilized 
rTMS of the DLPFC in major depression were included in the current meta-analysis. The sensi- 
tivity and moderator analyses also included data from N=40 studies (published in 1997-2008) 
from our previous meta-analysis. The effect size (Cohen's d) in each study was the standardized 
difference in mean depression scores (on Hamilton Depression Rating Scale, Beck Depression 
Inventory, Montgomery Asberg Depression Rating Scale) from baseline to final (after last 
session) in rTMS compared to sham groups. 

Results: According to a random-effects model with inverse-variance weights, depression 
scores were significantly reduced after rTMS compared to sham in studies published from 
2008-20 1 3 based on N=659 patients (overall mean weighted d=-QA2, 95% confidence interval: 
-0.66, -0.18, P=0.001). Combining studies from our past and current meta-analyses (published 
in 1997-2013; N=54) revealed that depression was significantly reduced after left-fast (>1 Hz), 
right-slow (^1 Hz), and bilateral (or sequential) rTMS of DLPFC compared to sham. Significant 
antidepressant properties of rTMS were observed in studies with patients who were treatment 
resistant, unipolar (or bipolar), non-psychotic, medication-free (or started on antidepressants 
concurrently with rTMS). According to univariate meta-regressions, depression scores were 
significantly lower in studies with more female patients and fewer stimuli per session. There 
was little evidence that publication bias occurred in the analysis. 

Conclusion: According to this study, the largest meta-analysis to date, short-term antidepres- 
sant properties of rTMS are independent of concurrent antidepressants and might depend on 
sex and the number of stimuli per session. 

Keywords: repetitive transcranial magnetic stimulation (rTMS), depression, sham-controlled, 
DLPFC, meta-analysis, systematic review 

Background 

A large volume of academic publications has been dedicated to the antidepressant 
properties of repetitive transcranial inagnetic stimulation (rTMS) in the treatment of 
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major depression. Our search of the MedHne and Psyclnfo 
databases identified N=963 sources (duplicates excluded) 
with terms "rTMS" and "depression" in their titles or sub- 
ject between (any date to September 2013). A vast majority 
of these sources are narrative literature reviews largely 
suggesting that approximately ten sessions of daily rTMS 
appear to be effective in acute cases of major depression (or 
major depressive episode) in the short-term (by comparing 
depression scores before the first versus after the last session 
of rTMS). However, the exact factors mediating the antide- 
pressant properties of rTMS are still not well understood. 
According to randomized-controlled trials (RCTs) and 
open-label studies conducted on mostly unipolar patients, 
rTMS was more effective in patients who were younger,''^ 
less treatment-resistant (in the current episode or with less 
prior treatment failures),'"'"^ with a shorter current episode,^-' 
and without a comorbid anxiety disorder.^ Furthermore, 
extension trials in patients who failed to respond during 
the short-term, double-blind phases of studies showed that 
the antidepressant response to rTMS was superior in female 
patients,^ was observed only after longer stimulation periods 
(such as 4 weeks or more),^'* and required an alteration in the 
stimulation site and frequency (from left-fast to right-slow 
rTMS of dorsolateral prefrontal cortex [DLPFC]).' 

Surprisingly, the antidepressant predictors of rTMS from 
the primary studies listed above have not been systematically 
confirmed in the relevant 17 quantitative meta-analyses 
(published 2001-2013) of the high-quality primary studies 
(sham-controlled randomized trials). '"^^ According to the 
meta-analyses to date, the short-term antidepressant prop- 
erties were most consistently observed in studies using the 
fast (> 1 Hz) rTMS of the left DLPFC.^^-^'' The slow (< 1 Hz) 
rTMS of the left or right DLPFC and bilateral or sequential 
designs were also effective at reducing depression severity 
in the short-term but were utilized in only very few primary 
studies. '^-"^-'''^'-^^ 

Similarly to the RCTs, the benefit of longer study designs 
(with ten or 15 rTMS sessions) in treating depression has 
already been noted in the earlier meta-analyses.' '^ '^ However, 
neither duration of study and other rTMS parameters (fre- 
quency of stimulation, motor threshold, stimuli/session, total 
stimuli) nor mean age of patients were associated with the 
effect sizes in meta-analyses. ^ ''''"•"■^''■^^ Furthermore, rTMS 
was effective in studies with medicated or medication-free 
patients, "•'''•"■^^'^' as well as in studies with medication- 
resistant patients. ''•■"'•"•^'' However, a better outcome was 
expected with less resistance.*" Finally, the antidepressant 
effect of rTMS was higher in studies with non-psychotic 



patients" but was similar in studies with unipolar versus 
bipolar patients. ^"-^^ 

One reason for such inconsistent findings is that most 
past meta-analyses included too few studies to reliably detect 
any differences in effect sizes based on study characteristics 
(clinical and/or rTMS parameters). Furthermore, unlike in 
meta-analyses, predictors of rTMS response were often iden- 
tified during different (open-label and/or follow-up) phases 
of primary studies. Finally, meta-analyses were computed 
based on group data compared to primary studies that had 
utilized individual patient data. 

In an attempt to improve the statistical power of the 
past analyses, we have conducted a meta-analysis on N=40 
sham-controlled studies selected from the past 13 meta- 
analyses'*''"''''"'"" published between 2001 and 20 10. A short- 
term antidepressant effect of the left-fast rTMS of DLPFC was 
univariately observed in studies with higher proportions of 
female patients not controlling for any other study character- 
istics (clinical and/or properties of rTMS). The antidepressant 
effect of the left-fast rTMS was also present in studies with 
patients who were medication-free, unipolar (or bipolar), 
treatment-resistant and without psychotic features. 

The current study had three main aims. Since our previ- 
ous meta-analysis included primary studies published up to 
2008, the first aim of the current study was to update our 
results by conducting a new meta-analysis of the short-term 
effects of rTMS in depression in studies published after 
2008 until August 2013. These "new" studies were located 
using a novel systematic literature search in contrast to 
the N=40 "old" studies in our previous meta-analysis that 
were selected from the past 13 meta-analyses published in 
200 1-20 1 0. Thus, the second aim of the current study was to 
compare the overall mean weighted effect sizes of the "old" 
studies with the "new" studies due to the different methods 
of searching for primary studies utilized in the two meta- 
analyses. Furthermore, our previous meta-analysis focused 
on the characteristics of studies that had utilized only the 
left-fast rTMS of DLPFC. Thus, the third aim of the cur- 
rent analysis was to find out if any patient characteristics or 
rTMS parameters would be associated with the short-term 
antidepressant properties of rTMS in all "new" and "old" 
sham-controlled studies published between September 1997 
and August 2013. The reason for combining all studies was 
to improve the statistical power of all statistical (moderator 
and subgroup) analyses. 

Based on our and other past meta-analyses, it was hypoth- 
esized that depression would be reduced following the active 
rTMS compared to sham in the "new" studies (those after 
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2008). It was expected that such an antidepressant effect 
would be higher in the "new" compared to the "old" studies 
if the quality of the more recent studies has improved due to 
more advanced stimulators and better established parameters 
of rTMS. Based on our past meta-analysis of the left-fast 
rTMS studies, it was expected that, when combining all stud- 
ies regardless of rTMS parameters, depression scores would 
be significantly reduced in studies with higher proportions 
of female patients. We also expected that, based on results 
from primary studies, the antidepressant properties of rTMS 
could be related to other patient characteristics and/or rTMS 
parameters if the statistical power of such comparisons were 
improved by adding the "new" studies to the "old" ones. 

Methods 

Systematic literature search 
and study selection 

The details of the systematic literature search are shown 
in Table 1 . A "control search" was first conducted in the 
Psyclnfo and Medline databases for N=40 studies pub- 
lished between 1995 and 2008 that were included in our 
previous meta-analysis. These N=40 studies were obtained 
from the past 13 meta-analyses (published 2001-2010) 
rather than from a systematic literature search (Table 
SI). Since all N=40 studies were located during the con- 
trol search, we concluded that these two databases were 
adequate for performing the current literature search for 
studies published in (any month of) 2008 until August 
2013. 

The results of the systematic literature search and the 
study selection procedure are summarized in the PRISMA 



flowchart (Figure 1).^^ Following the exclusion of irrelevant 
studies (based on titles and abstracts), N=50 primary studies 
were assessed in full-length (Figure 1 and Table S2). A total of 
N=18 out of 50 "new" studies (pubhshed2010-2013) located 
during our systematic search met the inclusion criteria for the 
current analysis (none of the studies published in 2009 met 
the inclusion criteria). Most studies were excluded because 
they were not sham-controlled or contained data published in 
other studies already included in the current analysis (other 
exclusion criteria are listed in Figure 1 ). The inclusion criteria 
for the current meta-analysis were; 

1 . sham-controlled parallel design; 

2. major depressive disorder or episode diagnosed 
according to the Diagnostic and Statistical Manual of 
Mental Disorders (DSM-IV) or International Statistical 
Classification of Diseases and Related Health Problems 
(ICD-10) criteria; 

3. depression severity assessed using any version of 
a standardized scale (Hamilton Depression Rating 
Scale (HAMD),^* Beck Depression Inventory (BDI)," 
and Montgomery Asberg Depression Rating Scale 
(MADRS^'*)); 

4. and active rTMS and sham administered at the same 
location of DLPFC (left, right, or bilateral). 

As explained in the Results section, four of the 1 8 "new" 
studies were identified as outliers and removed from all 
analyses. Thus, the current meta-analysis was performed on 
N=14 "new" studies. To improve the power of comparisons, 
the sensitivity analyses also included data from N=40 "old" 
studies (published in 1997-2008) from our previous meta- 
analysis (Table SI). 



Table 1 Details of the systematic search strategy (all searches were performed in English with 


no language restrictions) 


Search 


Search terms 


Databases (time frame) 


Search 1 


[Tl or SU ("transcranial magnetic stimulation" or "trans- cranial magnetic 


Psyclnfo and Medline 


N=82l 


stimulation" or "repetitive transcranial magnetic stimulation" or "repetitive 


(any date - 2008) 




trans- cranial magnetic stimulation" or TMS or rTMS)] AND [Tl or SU 






(depress* or dysthy* or MDD or cyclothym* or bipolar or "mani*- depress*")] 




Search 2 


[Tl or SU ("transcranial magnetic stimulation" or "trans- cranial magnetic 


EBSCO Psyclnfo and Medline 


N=S84» 


stimulation" or "repetitive transcranial magnetic stimulation" or "repetitive 


(2008 -August 20 13) 




trans- cranial magnetic stimulation" or TMS or rTMS)] AND [Tl or SU 






(depress* or dysthy* or MDD or cyclothym* or bipolar or "mani*- depress*")] 




Search 3 


[Title, Abstract, Keywords ('transcranial magnetic stimulation" or "trans- 


Cochrane Library (Category 


N=I28' 


cranial magnetic stimulation" or "repetitive transcranial magnetic stimulation" 


searched: Trials) 




or "repetitive trans- cranial magnetic stimulation" or TMS or rTMS)] AND 


(2008 -August 20 13) 




[Title, Abstract, Keywords (depress* or dysthy* or MDD or cyclothym* or 






bipolar or "mani*- depress*")] 





Notes: Search I was a "control search" to find out if the N-40 studies included in the past I 3 meta-analyses (published in 200 1— 20 1 0) could be located using two databases 
only. Since the search detected all these resources, Search 2 was conducted on Psyclnfo and Medline that appeared to have an adequate coverage of studies in this area. 
Search 3 of the Cochrane library did not identify any additional studies than Search 2. ^Duplicates excluded within search. 

Abbreviations: MDD, major depressive disorder; N, number of sources; rTMS, repetitive transcranial magnetic stimulation; SU, subject; Tl, title. 
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Identification 



Screening 



N=712 "new" records 
published in 2008- 
2013 from database 
searching 




N=40 "old" studies published in 
1997-2008 from our previous 
meta-analysis 








N=585 "new" records after duplicates removed 



N=585 "new" records 
(titles/abstracts) 
screened by VA/KKK 






N=535 excluded 







Eligibility 



Included 



N=50 "new" full-text 
articles assessed by 
VA/KKK 



N=18 "new" studies 
included in the 
quantitative meta- 
analysis (some 
analyses included 
N=40 "old" studies 
from our previous 
meta-analysis) 



N=32/50 (64%) articles 
excluded: 

• N=14: No sham group 

• N=10: No new data (data 
from studies already in 
analysis) 

• N=3: Inadequate data 
reported and lack of 
response from authors 

• N=1: Cross-over design 

• N=1: Case study 

• N=1: Depression secondary 
to pregnancy 

• N=1 : DLPFC not stimulated 

• N=1: rTMS after sleep 
deprivation 



Figure I Study assessment and exclusionary criteria. 

Abbreviations: DLPFC, dorsolateral prefrontal cortex; rTMS, repetitive transcranial magnetic stimulation. 



Data extraction 

Data were extracted from all N=18 studies by two authors 
(VA and KXK) independently and any inconsistencies were 
resolved (there were no major inconsistencies requiring addi- 
tional experts' opinion). The rTMS parameters are shown in 
Table 2, and the clinical and demographic characteristics of 
patients are shown in Table 3. 

Meta-analysis 

The mathematical approach used in the current 
meta-analysis is based on the method of Hedges et al.^' 
The mean depression scores at baseline as well as after 



the last session (final) in the sham and the active rTMS 
groups in each study are listed in Table 3. The mean (M) 
and standard deviations (SD) of depression scores in the 
sham and the active rTMS groups were computed for each 
group separately in SPSS version 21 (IBM Corporation, 
Armonk, NY, USA) using the difference score (depression 
score at baseline - final). These difference scores and their 
SD are shown in the last two columns of Table 3. Multiple 
independent subgroups of patients were combined into one 
active rTMS and one sham group per scale and per study 
to comply with the assumption of meta-analysis that each 
study should contribute only one effect size to the overall 
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analysis. These combined scores appear in rows labelled 
"all" in Table 3. 

The meta-analysis was computed using Comprehen- 
sive Meta-Analysis 2.0 (CMA; Biostat Inc., Englewood, 
NJ, USA). If one study utilized multiple depression scales 
then multiple effect sizes based on each scale were com- 
puted; these effect sizes were combined into one effect 
size per study using an arithmetic mean. The effect size 
used in the current analysis was the standardized mean 
difference (Cohen's d) between depression scores at base- 
line to final after sham compared to active rTMS. The 
interpretation criteria for the absolute size of Cohen's d 
are: £/=0.20-0.49 (small), t/=0.50-0.79 (moderate), and 
t/>0.80 (large).3° 

The current study utilized the random-effects model of 
meta-analysis^' with the inverse-variance weights (a sum 
of the within- and between-study variance)." The overall 
mean weighted effect size d of all studies was the sum of the 
product of all effect sizes and weights divided by the sum of 
all weights.^' Negative values of d indicated that rTMS had 
antidepressant properties compared to sham. 

The heterogeneity among the effect sizes was investigated 
using a Q statistic and an F index (F=\QQ%-x.{Q-df)IQ with 
df=k-\; i-number of studies).^' The P index shows the 
variability in effect sizes due to real (rather than chance) 
differences among studies and can be interpreted using the 
following criteria: 25% (little heterogeneity due to real dif- 
ferences among studies), 50% (moderate heterogeneity), and 
75%) (high heterogeneity).'^ 

Sensitivity and moderator analyses 

The stability of the overall mean weighted d over time was 
investigated as one study at a time was added to all previ- 
ous studies (cumulative analysis) and as one study at a time 
was removed from the overall analysis (one study removed 
analysis). The moderator analyses (subgroup analyses and 
univariate meta-regressions) were used to test the influence of 
systematic differences among studies (clinical characteristics 
or patients and rTMS parameters) on the overall mean 
weighted d. 

Publication bias analyses 

Publication bias was assessed using methods available 
in CMA (Biostat). The Rosenthal's Fail-Safe N" was 
computed to find out the number of studies (theoretically 
missing from the current analysis) required to lower the 
overall mean weighted d to zero in the current analysis. 



The Duval and Tweedie's Trim-and-Fill analysis'" was used 
to test if the so-called funnel plot of d versus standard 
error of the mean (SEM)/study'^ was symmetrical around 
the overall mean weighted d of all studies. Finally, the 
Begg and Mazumdar Rank Order Correlation (Kendall's 
tail b) between the standardized d versus SEM in each 
study,'' and the Egger's regression of 1/SEM (predic- 
tor) on the standardized d, " were used to test whether 
smaller studies differ systematically (significantly) from 
the larger studies. It was assumed that publication bias 
might be present if the Fail-Safe N is low, the funnel plot is 
asymmetrical, Begg and Mazumdar Correlation is statisti- 
cally significant, and/or the intercept of the regression line 
significantly deviates from zero, causing an asymmetry of 
the funnel plot.^' 

Results 

An inspection of the weighted effect sizes d for each of 
the N=58 studies revealed that N=4 studies (all "new" 
studies)'^""* were outliers in the current analysis. Specifi- 
cally, weighted <is in these studies were statistically sig- 
nificantly higher than the overall mean weighted d of all 
other studies (Figures SI and S2). As a consequence, the 
overall mean weighted d was inflated when these studies 
were included in the meta-analysis of the N=18 "new" 
studies ((i=-0.80) compared to when these studies were 
removed from the analysis (N=14 "new" studies: d=-QA2; 
Table 4). Thus, these N=4 "new" studies were excluded 
and all subsequent analyses were computed using N=14 
"new" studies. 

One of the four studies" might have been an outlier 
because of the following methodological differences between 
this and all other studies: the use of a clinical interview for 
HAMD scale (rather than HAMD scale alone), depression 
diagnosed using ICD-10 (rather than DSM-IV), inclusion 
of high proportions of patients with bipolar (30%)) and psy- 
chotic depression (68%; Table 3). The reasons why the other 
three studies might have been outliers are addressed in the 
Discussion section. 

The N= 14 "new" studies were conducted in nine countries 
(three each in the US and People's Republic of China; two 
in Spain; one each in France, Italy, Turkey, India, Canada, 
and Australia; Table 2). These studies were conducted on 
a total of N=659 patients in the active rTMS (N=340) and 
sham (N=3 19) groups. All patients had diagnoses of a major 
depressive episode and/or disorder according to DSM-IV 
or ICD-10 (Table 3). There was little evidence that the 
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Table 4 Results of the random-effects meta-analysis of the N=I4 "new" studies (2010-2013) and all N=54 studies: N=40 "old" studies 
(1997-2008) and N=I4 "new" studies 





N studies (%) 


d (9S% CI) 


p 

two-tailea 


V^/VCl d.11 1 1 ICdl 1 WClgllLCU U 








"old 40" studies (1997-2008) 


40 


-0.54 (-0.68, -0.41) 


<0.00l* 


"npw 1 ft" <:i-i iHipi on 1 f)— 9n i '^^ 


18 


fi on f \ \ fi. n AA\ 

— U.oU ^— 1 . 1 D, — U.T-'r ) 


nn i * 

*^U,UU 1 


"new 14" studies, outliers excluded 


14 


-0.42 (-0.66, -0.18) 


<0.00l* 


"old 40" versus "new 14" studies 






0.151 


"New" N= 1 4 studies 


90 1 0—90 1 ^ 






Total N patients 


rTMS/sham 


659 (340/319) 




Heterogeneity 




Q(df I3)=25.9; 
P=0.0 1 8*; P=50% 




Rosenthal's Fail-safe N for P>0.0S 


N— studies needed to reduce the 
overall mean weighted d to 0 


N=7I (N=5 studies missing 

tot" pvprv iHv inf"liiHfin in 

the current analysis) 




LJuval and 1 weedie s trirn-and-tii! analysis 


runnel plot symmetrical; 

N studies missing on either side 

of the overall mean weighted d 


1 es 
None 




Begg and Mazumdar rank order correlation 


T, P 

' two-tailed 


-0.09; P=0.66l 




Egger's regression 


Intercept; P 


-1.21; P=0.276 




Note: *P<0.05. 








Abbreviations: CI, confidence interval: d, standardized 


mean difference (effect size): df, degrees of freedom: rTMS, 


repetitive transcranial magnetic stimulation; 


r, Kendall's 



correlation coefficient tau b v^ith continuity correction. 



publication bias systematically affected the results of the 
current analysis (Table 4). 

There was a moderate antidepressant effect of rTMS 
because the change in mean depression scores from baseline to 
final was significantly higher after rTMS compared to sham in 
the N=14 "new" studies (c?=-0.42; Table 4 and Figure 2). The 
"new" studies did not show a superior antidepressant effect of 
rTMS compared to the "old" studies because the overall mean 
weighted d did not significantly differ between the two groups 
of studies ("old" d=-0.54 versus "new" J=-0.42; P=0.151; 
Table 4). 

The 50% heterogeneity in the effect sizes among the N=14 
"new" studies (Table 4) was probably due to methodologi- 
cal differences among these studies in terms of depression 
scales used (HAMD in all N=14 studies, MADRS in N=6 
studies, and BDI in N=3 studies), clinical characteristics 
of patients (Table 3), and different parameters of rTMS 
(Table 2). Some of the N=14 "new" studies did not report 
the above characteristics (Tables 2 and 3) and thus the per- 
centage scores shown below are computed based on studies 
with valid responses only. 

The N=40 "old" studies and the N= 1 4 "new" studies were 
conducted on mostly middle-aged patients (aged 40 years old 
and above). Similarly to the "old" studies, 50% of the "new" 
studies utilized 60% or more female patients per study. 



Like the "old" studies, most of the "new" studies were 
conducted on patients with (Table 3): 

• treatment-resistance, defined as a failure to respond 
to or tolerate s2 antidepressant trials (N=8/9, 89% of 
studies), 

• non-psychotic depression (N=8/9, 89% of studies), and 

• concurrent antidepressant treatment (N=ll/14, 79% of 
studies). Of these N=l 1 studies, most included patients 
on stable doses of antidepressants (N= 10/11 studies). 
Antidepressants were started on day 1 of the study in 
only N= 1/11 studies. 

In contrast to the majority of "old" studies that had 
included any proportions of bipolar patients, the "new" 
studies were mostly conducted on patients with unipolar 
depression (N=8/12, 67% of studies). 

The current results also suggest that similar properties of 
rTMS have been used over the last 16 years of research on 
depression (in studies published in 1997-2013). Similarly to 
the "old" studies, the most commonly utilized rTMS param- 
eters among the "new" studies were: 10 Hz frequency of 
stimulation (N=5/14, 36% of studies), 1 10% motor threshold 
stimulation (N=4/13, 41% of studies), 800 or 1,600 stimuli/ 
session (N=4/10, 40% of studies), 20 trains/session (N=6/13, 
46% of studies), a 70 mm stimulating coil diameter (N=4/5, 
80% of studies), and a figure-of-eight shape of the stimulating 
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statistics for each study 



Sample size 



Std diff in means and 95% CI 







Std diff 


Lower 


Upper 












in means 


limit 


limit 


P-value 


Sham 


rTMS 


Oif^nrnf^ fit al^^ 9010 

OcUiyc al lU 


womDi neo 


-Oil 


-0 40 


0 18 


0 456 


94 


88 


Paillere Martinet et a!^" 2010 


Combined 


-0.40 


-1.10 


0.31 


0.269 


14 


18 


Pallanti et al^'' 201 0 


HAMD 


-1.33 


-2.01 


-0.64 


0.000 


20 


20 


Triggs et al^^ 2010 


Combi ned 


-0.19 


-0.81 


0.43 


0.548 


14 


34 


Am lirrfi f^i 901 1 

/Ayiillic CL al 1 1 


HAMD 


-0 18 


-0 86 


0 50 


0 605 


1 5 


1 9 


Ufa fat al^'t Oni -1 
rnc cL al 1 1 


nMIVILJ 


-0.44 


-0.88 


0.01 


0.054 


43 


37 


l—ingeswaran et al^^ 201 1 


Combi ned 


-0.09 


-0.93 


0.75 


0.836 


14 


g 


baKim et ai zui^; 


Combined 


—0.93 


—1 .79 


—0.07 


0.035 


1 2 


1 1 


Blumberger et al" 2012 


HAMD 


0.22 


-0.43 


0.87 


0.505 


18 


19 


Fitzgerald^^2012 


Combined 


-0.39 


-1.03 


0.24 


0.227 


18 


21 


Huang et a|5^2012 


Combined 


-0.95 


-1.51 


-0.40 


0.001 


28 


28 


Chen etaP 2013 


Combined 


0.38 


-0.55 


1.30 


0.427 


10 


10 


Hernandez-Ribas et al^^ 


HAMD 


-0.89 


-1.79 


0.01 


0.052 


11 


10 


Speer etaF2013 


HAMD 


-0.90 


-1.78 


-0.01 


0.047 


8 


16 






-0.42 


-0.66 


-0.18 


0.001 







rTMS effective Sham effective 



Funnel plot of standard error by std diff in means 




Std diff in means 



Figure 2 Random-effects meta-analysis of N— 1 4 "new" studies (20 1 0-20 1 3) comparing the change in mean depression scores on HAMD, BDI, and/or MADRS (baseline - 
final), after rTMS versus sham. 

Notes: In the studies by Triggs et al" and Speer et al," rTMS was administered using different properties into different subgroups of patients in a study and the depression 
scores for such subgroups were combined. The mean number of patients per group was used in the final calculations if patients dropped out throughout the study between 
baseline and final sessions. The forest plot (top) shows the weighted effect size d (box) and its 95% CI (vertical line through the box) for each study in the analysis ("combined" 
indicates that more than one depression scale was used in a study and the effect sizes according to the multiple scales were combined). The diamond depicts the overall 
mean weighted d of all studies and its 95% CI (width of the diamond). The mean depression scores (baseline - final) were significantly reduced after rTMS compared to sham 
(overall mean weighted d--0.41, 95% CI: -0.64, -0.18). The funnel plot (bottom) shows the effect sizes d plotted versus SEM for each study in the analysis. The plot was 
symmetrical around the overall mean weighted d suggesting that publication bias had little effect on the results of the current meta-analysis. 

Abbreviations: BDI, Beck Depression Inventory: CI, confidence interval: HAMD, Hamilton Depression Rating Scale; MADRS, Montgomery Asberg Depression Rating Scale: 
rTMS, repetitive transcranial magnetic stimulation: SEM, standard error of mean: Std diff, standardized mean difference d. 



coil (N=13/13 studies). In contrast to the 10 rTMS sessions 
most commonly used among the "old" studies, the "new" 
studies most often utilized longer paradigms of 15 sessions 
(N=6/14, 43% of studies). 

Interestingly, even though the new "sham-coils" improve 
the blinding of studies and prevent stimulation of the brain due 
to an inbuilt magnetic shield these coils have not been com- 
monly used in the "new" studies yet (possibly due to a high 
cost of replacing the older stimulators with the newest ones). 
In fact, similarly to the "old" studies, the most common sham 
practice among the "new" studies was to tilt the active coil by 
a 90° angle from the scalp (N=6/14, 43% of studies). 

Finally, just like the "old" studies, the left-fast rTMS 
(> 1 Hz) of the DLPFC was the most common combination 
of the location-frequency of rTMS among the "new" stud- 
ies (N=10/14, 71% of studies), followed by the bilateral or 
sequential rTMS (N=3/14, 21% of studies), and the left-slow 



(< 1 Hz), right-slow, and right-fast rTMS that were utilized 
in one study (7%) each. 

Since there were too few "new" studies, we have com- 
bined all "old" and "new" studies (total N=54) to conduct 
moderator and subgroup analyses. The results of these 
analyses are shown in Table 5. 

The overall mean weighted d=-0.5l in all N=54 studies 
suggests that rTMS was only moderately better than sham at 
reducing depression scores over the last 1 6 years in a total of 
2,242 patients in 17 countries (Table 5; Figure S3; Table S3). 
Even though only moderate, this overall effect size became 
stable around -0.50 to -0.54 over the last 8 years as studies 
were removed from the analysis one at a time (Figure S4), or 
cumulatively added to the analysis (Figure S5). 

Univariate comparisons of subgroups of studies based 
on common study properties showed that no one specific 
characteristic was superior in terms of producing higher 
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Table 5 Random-effects subgroup analyses and meta-regressions of 
compared to sham in N=54 sham-controlled studies published in 1997 



the change in depression scores (baseline 
- August 20 1 3 



final) after rTMS 



Studies 



N studies (%)' 



d (95% C\) 



All studies (1997-2013) 
Subgroups 

Location-frequency*" 

Left-slow (< I Hz) 

Left-fast (> I Hz) 

Right-slow 

Right-fast 

Bilateral or sequential (left then right) 
Treatment resistance 

Yes (all failed 2:2 AD trials) 

No (all failed 0-1 AD trials) 

Yes versus no' 
Concurrent medication 

YES (any % of patients) 
Stable dose 
Started on day I 

NO (all patients) 

YES versus NO' 
Bipolar depression 

YES (any % of patients) 

NO (all patients) 

YES versus NO' 
Psychotic depression 

YES (any % of patients) 

NO (all patients) 

YES versus NO' 
Coil-type 

F8 

Circular 

F8 versus circular' 
Coil angle sham 
0° (inactive coil) 
0° (sham coil) 
45° 
90° 

0° (sham coil) versus 45° 
0° (sham coil) versus 90° 



54 
54 

4 (7%) 
43 (80%) 
6(1 1%) 

I (2%) 

7(13%) 

27 

21 (78%) 
6 (22%) 

54 

42 (78%) 
35 (83%) 
7(17%) 
12 (22%) 

42 

23 (55%) 
19 (45%) 

28 

5 (18%) 
23 (82%) 

51 

47 (92%) 

4 (8%) 

54 

5 (9%) 
10 (19%) 
18 (33%) 
21 (39%) 



-0.5 1 (-0.63, -0.39) 



-0.61 
-0.49 
-I.OI 
0.03 ( 
-0.55 

-0.52 
-0.80 



-0.51 
-0.51 
-0.50 
-0.56 



-0.44 
-0.54 



-0.51 
-0.58 



-0.52 
-0.62 



-0.36 
-0.63 
-0.40 
-0.56 



1.21,-0.01) 
0.63, -0.34) 
1.61,-0.42) 
0.86, 0.92) 
-0.82, -0.29) 

-0.70, -0.35) 
-1.02, -0.50) 



-0.63, -0.38) 
-0.65, -0.36) 
-0.77, -0.23) 
-0.84, -0.28) 



-0.60, -0.28) 
-0.72, -0.34) 



-1.14, 0.13) 
-0.77, -0.40) 



-0.65, -0.38) 
-1.05,-0.19) 



-0.64, -0.07) 
-0.92, -0.34) 
-0.57, -0.22) 
-0.77, -0.35) 



N studies 



<0.00l* 



0.046* 

<0.00l* 

0.001* 

0.948 

<0.00l* 

<0.00l* 
<0.00l* 
0.108 

<0.00l* 
<0.00l* 
<0.00l* 
<0.00l* 
0.229 

<0.00l* 
<0.00l* 
0.921 

0.1 17 

<0.00l* 
0.745 

<0.00l* 

0.005* 

0.561 

0.015* 

<0.00l* 

<0.00l* 

<0.00l* 

0.757 

0.150 



Meta-regression predictors'* 
% female patients 
Stimuli/session 
Trains/session' 



53 
33 
48 



0.046 
0.043 
0.071 



0.022 
0.004 
0.044 



0.52 
0.91 
0.38 



-0.01 

0.0002 

0.007 



0.002* 

<0.00l* 

0.013* 



Notes: Total patients in N-54 studies totaled N-2,242 (rTMS N-I,I84, sham N-1,058). The percent values are reported based on the number of studies that reported 
a particular characteristic: ^'effect sizes in subgroups based on location-frequency of rTMS were not compared statistically because some studies used multiple active rTMS 
groups but the same sham groups, and thus the subgroups were not independent; 'subgroups were compared using the mixed-effect model; random-effects model was used 
to compute the overall mean weighted d in each subgroup and overall mean weighted d of subgroups were compared using the fixed-effect model because the number of 
subgroups was fixed; ^proportion of the between-study variance in weighted d explained by the predictor was computed as R^-l — {^^^oje/^ioid)' where T^^^^^, is the between- 
study variance in the weighted d (outcome) unexplained by the regression model containing the predictor and T^^^^^^ is the within- and between-study variance;^' 'following the 
Bonferroni correction for multiple regressions {new significance threshold of 0.05/7—0.007), the regression of trains/session on weighted d became non-significant. *P<0.05. 
Abbreviations: AD, antidepressant; CI, confidence interval; rTMS, repetitive transcranial magnetic stimulation. 



antidepressant effects. For example, except for the right- 
fast rTMS of the DLPFC that was utilized in only one 
study, all combinations of the location-frequency of rTMS 
were effective at significantly reducing depression scores 
compared to sham (Table 5). There was a non-significant 



trend toward higher overall mean weighted effect sizes in 
studies with non-treatment resistant patients compared to 
those with treatment-resistant patients (Table 5). rTMS was 
also similarly effective in studies with unipolar depression 
compared to studies with generally low proportions of 
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bipolar patients (Table 5). However, it remains unclear if 
the antidepressant properties of rTMS extend to psychotic 
depression because there were too few studies with low pro- 
portions of patients with psychotic depression in the current 
analysis (N=5; Table 5). 

According to the current results, the antidepressant 
effect of rTMS was probably not secondary to concurrent 
antidepressants. Table 5 shows that depression was reduced 
after rTMS compared to sham in studies with patients 
who were medication-free or started on antidepressants 
concurrently with rTMS on day 1 of a study. Overall mean 
weighted effect sizes did not differ statistically between 
studies with medicated versus medication-free patients 
(Table 5). 

The current study also shows that depression was reduced 
after rTMS compared to sham using both stimulating coil 
shapes (figure-of-eight or circular) and various sham designs. 
Although studies using sham coils produced the highest over- 
all mean weighted effect size (t/=-0. 63; Table 5) compared to 
all other sham designs, more than N=10 studies using sham 
coils are needed to statistically confirm the superiority of this 
blinding method. Compared to sham coils, tilting of active 
coils at 45° or 90° from the scalp was the most commonly 
used sham practice in studies published until August 2013 
(Table 5). The most commonly used stimulators were the 
MagStim (UK; N=28/54, 52% of studies) and the MagPro 
(USA; N=10/54, 18% of studies) models. 

Finally, seven univariate meta-regressions were conducted 
to find out if any demographic characteristics of patients 
(mean age/study and proportion of female patients/study) or 
rTMS parameters (frequency of stimulation, motor threshold, 
total number of sessions, stimuli/session, trains/session) 
could predict a change in the effect sizes weighted according 
to the random-effects model. Two regressions showed that the 
antidepressant effect of rTMS was superior in studies with 
more female patients and fewer stimuli per session (Table 5; 
Figure 3). However, a change in weighted ds could not be 
significantly predicted using the mean age of all patients per 
study, frequency of stimulation, motor threshold, or a total 
number of sessions as predictors. Predictors "stimuli/session" 
and "% female patients" explained 91% and 52% of the 
between-study variability in effect sizes, respectively. 

Discussion 

The current meta-analysis conducted on N=54 studies pub- 
lished in 1997-2013 showed that rTMS has a short-term 
antidepressant effect that is superior in studies with more 
female patients and fewer stimuli per session (Table 5 and 



Figure 3). The most commonly used characteristics of 
patients and rTMS parameters in the sham-controlled studies 
over the last 16 years (1997-2013) included in the current 
analysis were: at least 60% female patients/study in half of 
all studies (range: 22%]-92%); mostly middle-aged or older 
patients (range of the mean age of all patients per study: 
39-62 years); fast (>1 Hz) stimulation of the left DLPFC; 
frequency of 10 Hz; 1 10% motor threshold; figure-of-eight 
coil with 70 mm diameter; 1,600 stimuli/session; 20 trains/ 
session; and 10 sessions/study (although the more recent 
studies published after 2008 most commonly used longer 
protocols of 15 sessions/study). The most commonly used 
sham strategy was tilting of the active coil at the 90° angle 
from scalp (Table 5). Most of the studies to date included 
proportions of patients with treatment resistance, on concur- 
rent antidepressants (particularly at the stable dose), with 
bipolar and non-psychotic depression (Table 5). 

It is likely that sex plays a role in the short-term 
response to the left-fast rTMS because 80% of all studies 
in the current meta-analysis utilized this combination of 
rTMS parameters (Table 5). This result confirms another 
univariate finding that out of patients who failed to respond 
to a 4-week, double-blind phase of rTMS in a large RCT,"^ 
only females showed a superior response to rTMS during the 
extension (open-label) phase of the study.^ However, such a 
response to treatment was probably related to a combination 
of factors rather than sex alone. This is because the patients 
in the RCT were unipolar, non-psychotic, medication-free, 
moderately-severely treatment-resistant, and required more 
than four weeks of treatment to respond to rTMS.^ Another 
open-label study also showed that the improved response to 
rTMS in females depended on younger (premenopausal) age 
and the ovarian hormonal levels.''^ Such a result is not surpris- 
ing because depression has a strong hormonal component.'''' ''^ 
Thus, it can be speculated that the superior effect of sex in our 
meta-analysis was also due to yoimger age of female patients, 
and other factors (such as less severe unipolar depression 
and/or medication-free status). However, such confounding 
factors can be investigated only to a limited extent in meta- 
analyses that are computed on data. It should be possible to 
compute multiple meta-regressions on group data as new 
studies using rTMS become available in the future, providing 
that these studies report the characteristics of their patients 
and/or rTMS properties used. 

The second important finding in the current meta-analysis 
is that the short-term antidepressant properties of rTMS were 
observed in studies using fewer stimuli per session accord- 
ing to a univariate meta-regression (Table 5 and Figure 3). 
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Other meta-analyses showed that the efficacy of rTMS was 
not associated with the number of stimuli/session nor total 
stimuli. The initial negative correlation between the 
effect sizes and the total number of stimuli in one meta- 
analysis^" was attributed to one large RCT only.''^ Similarly, 
this and another large RCT''^ "'* contributed to the additional 
significant univariate meta-regression of the total number 
of stimuli (computed as "stimuli per session" x "total ses- 
sions" per study based on the data shown in Table 2) on 
weighted d in the current analysis (Figure S6). However, the 
meta-regression of stimuli/session on weighted d remained 
imchanged after the removal of the same two RCTs*^ '"' in the 
current analysis (Figure S7). In general, the findings fi-omthe 
two largest RCTs'*^ '*'' suggest that the initially non-responsive 
patients appear to indeed require more than 15 rTMS ses- 
sions with a high number (3,000) of stimuli/session to show 
a response to rTMS. However, the results of the current 
meta-analysis suggest that the short-term response to rTMS 
(during the double-blind phase of the study) might require 
fewer stimuli per session. Again, based on the current results, 
it can only be speculated that particularly the less treatment- 
resistant female patients require fewer stimuli per session in 
the short (up to 10 sessions), left-fast rTMS paradigms to 
demonstrate the antidepressant response to rTMS. 

Interestingly, some evidence in support of the specula- 
tion above can be found in three of the four studies clas- 
sified as outliers in the current analysis. Specifically, a 
large antidepressant effect of rTMS was observed in the 
total of N=86 patients (63%-65% males) on concurrent 
antidepressants, with unipolar, non-psychotic depression, 
and with moderate-severe treatment-resistance using long 
(20 sessions), lefl:-fast paradigms with a high number (3,000) 
of stimuli/session.^* '"''" Therefore, in contrast to female 
patients, male patients with more severe major depression 
might require longer, left-fast paradigms with more stimuli per 
session to show an antidepressant response to rTMS during 
the double-blind phases of studies. Furthermore, such a sex- 
dependent effect could also explain why the two largest RCTs 
to date, with similar study characteristics to those in the three 
outlier studies, have demonstrated only small antidepressant 
effects of rTMS in mostly female (medication- free) patients: 
d=-03f) in N=301 patients''^ (Figure SI) and d=-0.\\ in 
N=190 patients'"' (Figure 2). Such small effect sizes might 
have resulted from a high quality of blinding (with sham 
coils) in the RCTs compared to the large effect sizes using 
tilted active coils in the outlier studies. Thus, future primary 
studies should investigate the effects of rTMS separately in 
both sexes controlling for severity of treatment resistance 



and the number of stimuli/session in the left-fast rTMS 
paradigms. 

The significant reduction in depression scores after rTMS 
associated with fewer stimuli per session raises the question 
of whether or not the brain can be "overstimulated" during 
rTMS leading to a reduction in the antidepressant properties 
of this method. One mechanism of such an "overstimulation" 
could be related to the firing properties of neurons. In general, 
voltage-gated sodium charmels are key players in membrane 
excitation and the production of action potentials. The clas- 
sical model of sodium channel gating described by Hodgkin 
and Huxley suggests that a voltage-gated mechanism medi- 
ates the activation (opening of sodium channels) as well as 
inactivation following a refractory period during which no 
excitation can occur.'" A strong depolarization involving 
a large number of neurons can inactivate sodium channels 
and thereby prevent further excitation for a prolonged period 
of time. Since neurons need to recover from firing before 
being able to produce new action potentials, stimulating the 
brain with too many stimuli might lead to a neural saturation 
(inability of most neurons to produce new action potentials) 
and consequently a reduction in antidepressant properties of 
rTMS. In practical terms, using shorter sessions with fewer 
stimuli could be less costly as well as time consuming for 
patients and administrators. 

Univariately, the mean age of all patients (rather than 
age of individual patients) was not related to a better anti- 
depressant outcome of rTMS in the current and past meta- 
analyses'"-^" possibly due to the use of group data. Similarly, 
effect sizes were unrelated to severity of treatment-resistance 
or presence versus absence of treatment-resistant patients in 
other meta-analyses. '''•"'■"•^'' In our analysis there was only a 
trend toward higher antidepressant effect of rTMS in (very 
few) studies with non-treatment-resistant patients compared 
to studies with treatment-resistant patients (Table 5). The 
reason for this result could be that we have not controlled 
for severity of illness in our meta-analysis. Furthermore, 
our classification of studies into subgroups was also not 
optimal for unipolar versus bipolar and psychotic versus 
non-psychotic depression. In general, very few patients 
had such diagnoses per study (Table 3). Thus, there was 
no difference in effect sizes between groups of studies 
with all unipolar and all non-psychotic patients compared 
to studies with mostly unipolar and mostly non-psychotic 
patients, respectively (Table 5). Other meta-analyses have 
also found no differences in effect sizes or response and 
remission rates between studies with unipolar compared to 
bipolar patients.''''^' While effect sizes were higher in studies 
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Figure 3 Univariate random-effects meta-regressions of various study characteristics used as predictors (proportion of female patients/study and stimuli/session) on the 
weighted effect sizes d (the outcome) in studies published in 1997-2013. 

Notes: The figures are scatterplots of the outcome (weighted d/study; Y-axes) versus predictors (X-axes): (A) proportion of female patients/study and (B) stimuli/session. 
Depression scores (baseline — final) were significandy reduced after rTMS compared to sham in studies with more female patients and less stimuli/session. 
Abbreviations: rTMS, repetitive transcranial magnetic stimulation; Std diff, standardized mean difference d. 
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with non-psychotic versus psychotic patients in one meta- 
analysis," such result was not confirmed by another meta- 
analysis,''* possibly due to including studies with unknown 
psychosis status in the non-psychotic group. 

The short-term antidepressant effect of rTMS was not sec- 
ondary to concurrent antidepressants in the current analysis. 
Specifically, it was observed in studies with patients who were 
all medication-free or who started on antidepressants concur- 
rently with rTMS (Table 5). However, the current analysis 
did not confirm the finding from other meta-analyses that 
the short-term response to rTMS was higher^^ or tended to 
be higher" in studies using rTMS as monotherapy versus an 
add-on therapy. The reason for this result might be that we 



have not controlled for proportions of patients on concurrent 
antidepressants in the add-on studies. 

According to our results depression severity was reduced 
after rTMS compared to sham in studies using different 
rTMS properties, such as different stimulating coils (figure- 
of-eight or circular), different sham paradigms, and different 
combinations of the location (right or left) and frequency 
(slow or fast) of DLPFC stimulation (Table 5). Although 
blinding is facilitated by the use of shielded sham coils that 
resemble active coils visually and produce similar auditory 
effects while not stimulating the brain,''* these coils were 
still not commonly used compared to tilting of active coils 
(possibly due to high costs of replacing the older with newer 
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equipment). Furthermore, except for five studies, all other 
studies used the "5 cm rule" to define the position of the 
DLPFC (Table 2) in the current analysis. Even though it is 
so frequently used, the "5 cm rule" is less accurate than the 
10-20 EEG (electroencephalogram) system and the magnetic 
resonance image (MRI)-guided neuronavigation.'" There- 
fore, the antidepressant properties of rTMS could be further 
strengthened by the use of either MRI or the F5 location of 
the EEG system.'*' 

The moderate effect sizes in the current and most other 
meta-analyses on this topic could be related to statistical 
methods of computing the effect sizes and performing a 
meta-analysis. The magnitude of effect sizes in individual 
studies might largely depend on the blinding quality of stud- 
ies and the computation of effect sizes. Therefore, assuming 
that baseline depression scores are similar, the well-blinded 
studies could have smaller standardized differences in mean 
depression scores between rTMS and sham groups. This is 
because patients could respond to rTMS (due to its anti- 
depressant properties) and to sham (due to placebo effect) 
in well-blinded studies. In contrast, the poorly-controlled 
studies in which patients and/or administrators guess the 
treatment allocation could show larger effect sizes. This 
time patients could respond to rTMS (guessing that they 
receive the real treatment) but not to sham (guessing that 
they receive the inactive treatment). This scenario is pos- 
sible because patients can easily inform themselves about 
the method from the social media, for instance by watching 
industry-sponsored videos on YouTube directed toward the 
general population. 

Furthermore, the past meta-analyses in this field used two 
general approaches to computing the effect sizes: the mean 
depression scores were compared before versus after rTMS 
or sham (in some studies controlling for baseline depression 
scores) using standardized mean differences (Cohen's d or 
Hedges' g), or proportions of patients who remitted after 
rTMS or sham versus baseline were compared using the odds 
ratios. None of these approaches of computing effect sizes 
is ideal from the statistical point of view. Mean scores could 
be skewed in studies with small sample sizes (<30 patients). 
Thus, the "true average" depression score would be shown 
more accurately using the mode (most commonly occurring 
depression score in the rTMS versus sham groups) or the 
median rather than the skewed mean. On the other hand, 
the odds ratios are also problematic because they rely on 
classifying patients into two groups (non-depressed versus 
depressed) based on a subjective cut-off on a depression 
scale. As a result, a patient with a score of 17 on a specific 



scale might be classified as "non-depressed" while another 
one with a score of 18 might already fall into the "depressed 
group" regardless of such a low difference between their 
scores. Furthermore, the odds ratios computed from each 
primary study would need to be based on the same cut-off 
for presence/absence of depression to reliably combine the 
results of these studies in meta-analysis. Therefore, the "all 
or nothing" classification of patients according to the odds 
ratio might be reliable and valid when large differences 
among patients' scores occur and thus the group member- 
ship can be reliably justified. Despite the limitations above, 
the advantage of mean scores is that they show the severity 
of depression without needing a specific cut-off to classify 
patients into groups. 

Furthermore, different methods of weighing of studies 
and meta-analysis were utilized in the meta-analyses to date. 
The weights in the current meta-analysis were computed 
based on variance. Therefore, studies with higher variability 
of scores (due to rTMS being effective at reducing depression 
in only some, but not all, patients) had lower weights and 
thus lower influence on the overall mean weighted effect size. 
However, other methods of meta-analysis advocate the use 
of other weighing methods, such as the sample size, and cor- 
recting the effect sizes in individual studies for study-related 
artifacts before conducing any meta-analysis.^' Despite all 
the statistical differences, most of the past 17 meta-analyses 
(published 2001-2013) and the current meta-analysis report 
similar (moderate) weighted effect sizes (standardized mean 
differences or odds ratios). Therefore, it is likely that the 
true effect of rTMS is either only moderate or indeed even 
higher in clinical practice considering the statistical limita- 
tions described above. 

Although no strong evidence for publication bias was 
detected in the current study (Table 4), the sources included 
in our analysis were biased toward studies published in peer- 
reviewed journals and written in the English language. Such 
an apparent bias in selection of sources was not related to 
the systematic search strategy (Table 1). Our search was 
conducted for any type of resource (published or not pub- 
lished) in any language (because the authors of this study are 
multilingual speakers of six different languages). Instead, it 
appears that the majority of sources on this topic on Psyclnfo 
and Medline were indeed published in peer-reviewed jour- 
nals and written in English (or at least included a title and 
an abstract written in English). Therefore, such a linguistic 
bias is probably related to the fact that Psyclnfo and Med- 
line mostly store published sources and these sources are 
most often written in English. However, the results of our 
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analysis are generalizable beyond the English-speaking 
world because the N=54 studies were conducted in a total 
of 17 countries around the world (for the list see Table S3). 
According to the File Drawer Problem" studies with statisti- 
cally significant results are more likely to be published while 
those with non-significant results remain in "file drawers" 
and are never published. Even though such so-called "gray 
(unpublished) literature" was not included, only 50% of all 
studies in the current meta-analysis reached the traditional 
significance level (Figure S3). Furthermore, no checklist 
for the quality of studies was used in the current analysis. 
Instead, the quality of studies was assessed indirectly by 
weighing the effect sizes based on variability of scores 
within and between studies. It was assumed that studies 
with low variability of scores were of higher quality and 
thus contributed more weight to the overall mean weighted 
effect sizes and vice-versa. 

Conclusion 

In conclusion, the results of the current study suggest that 
the short-term antidepressant properties of rTMS are not 
secondary to conciu^rent antidepressants and might depend 
on sex and the number of stimuli per session. Depending 
on degree of treatment-resistance and age, male and female 
patients with unipolar depression might require paradigms 
with different properties (number of stimuli per session, 
total number of sessions, left-fast or other combinations 
of location-frequency of stimulation) to show comparable 
antidepressant effects. While the clinical efficacy of the fast 
rTMS of the left DLPFC seems to be widely accepted, the 
right-slow and bilateral or sequential paradigms appear to be 
promising alternatives in the short-term treatment of acute 
major depression. 
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Table SI A list of N=40 "old" studies on the association 
between rTMS from our previous meta-analysis (unpublished). 
These studies were included in sensitivity and moderator analyses 
in the current study 

Author, year, country 

George etal,' 1997, USA 

Avery etal,^ 1999, USA 

Kimbrell etal,' 1999, USA 

Klein etal," 1999, Israel 

Loo etal,' 1999, Australia 

Padberg et al,' 1 999, Germany 

Berman et al,' 2000, USA 

Eschweiler et al,° 2000, Germany 

George et al,' 2000, USA 

Garcia-Toro et al,'° 2001, Spain 

Garcia-Toro etal," 2001, Spain 

Manesetal,'UOOI,USA 

Boutros," 2002, USA 

Padberg et al,'"* 2002, Germany 

Fitzgerald et al,'' 2003, Australia 

Hoppner et al," 2003, Germany 

Loo et al," 2003 Australia 

Nahas et al,'» 2003, USA 

Buchholtz et al," 2004, Denmark 

Hausmann et al,™ 2004, Austria 

Holtzheimer et al,^' 2004, USA 

Kauffmann et al,'' 2004, USA 

Koerselman et al,'' 2004, the Netherlands 

Mosimann et al,'"" 2004, Switzerland 

Poulet et al," 2004, France 

Rossini et al," 2005, Italy 

Rumi et al," 2005, Brazil 

Su et al,'" 2005, Taiwan 

Avery et al," 2006, USA 

Fitzgerald et al,'" 2006, Australia 

Garcia-Toro et al," 2006, Spain 

Januel et al,'' 2006, France 

Anderson et al," 2007, UK 

Bortolomasi et al," 2007, Italy 

Herwig et al,'' 2007, Germany/Austria 

Loo et al," 2007, Australia 

O'Reardon et al," 2007, USA, Australia, Canada 

Stern et al,'» 2007, USA 

Bredau et al," 2008, Denmark 

Mogg et al,« 2008, UK 

Abbreviation: rTMS, repetitive transcranial magnetic stimulation. 



Table S2 A list of N=50 studies on the association between 
rTMS and depression assessed in full-length and reasons for 
exclusion from the current meta-analysis 



Author, year 



Included/reason for exclusion 



Aguirre et al,"' 201 I 

Avery et al,"' 2007 

Baeken et al,"' 20 1 0 
Bakimetal,"" 2012 
Bares et al,"' 2009 



Blumbergeretal,"'20l2 
Brakemeier et al,"' 2007 
Brakemeier et al,"° 2008 
Chen etal,"' 20 1 3 
Cohen et al,'" 2009 
Dell'Osso et al," 2009 
Fitzgerald et al," 2012 
Furtado et al," 2012 
Galletly et al,'" 2012 
George et al," 2010 
Hadleyet al," 2011 
He et al," 201 I 

Herbsman et al,'= 2009 

Hernandez-Ribas et al," 2013 
Herwig et al,'" 20 1 0 

Hoppner et al," 2010 

Hoy et al,"20l2 
Huang et al," 2008 
Huang et al," 2012 
Jacob et al," 2008 

Kozel et al,"20l I 

Kreuzer et al," 2012 
Lingeswaran et al,'" 20 1 I 



Included (additional data provided 
by authors) 

No new data (data from 

Avery et al," 2006) 

Cross-over design, only one session 

Included 

No sham (rTMS and placebo 

medication versus sham and 

venlafaxine) 

Included 

No sham 

No sham 

Included 

No sham 

No sham 

Included 

No sham 

No sham 

Included 

No sham 

Included (additional data provided 
by authors) 

No new data (data from 
Avery et al," 2006) 
Included 

No new data (data from 
Herwig et al," 2007) 
No new data (data from 
Herwig et al," 2007) 
No sham 
No sham 
Included 

Inadequate data reported 
(SD values missing on Figure I) 
No new data (data from 
O'Reardon et al," 2007) 
rTMS after sleep deprivation 
Included 
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Table S2 (Continued) 



Author, year 



Included/reason for exclusion 



Lisanby et al," 2009 

Myczkowski et al,™ 20 1 2 

Nongpiur et al," 201 I 

Paillere Martinet et al," 20 1 0 

Pallanti et al," 2010 

Peng et al,"'20l2 
Rayet al,"20ll 
Rosenquist et al," 
2013 

Schrijvers et al," 2012 

Schutter et al,™ 2009 

Schutter et al," 2010 

Simpson et al,™ 2009 

Spampinato et al,°' 
2013 

Speer et al,»" 2009 

Speer et al,"' 
2013 

Tamas et al,'" 2007 

Triggs et al,=^ 2010 
Trojaketal,'"20l I 
Ullrich et al,='20l2 

Zarkowski et al,"" 
2009 

Zheng et al,'"20IO 
Zheng et al,"'20IO 



No new data (data from 
O'Reardon et al," 2007) 
Depression secondary to birth 
(postpartum depression) 
No sham (primed all conditions 
with right I Hz stimulation) 
Included (additional data provided 
by authors) 

Included (week 3 data extrapolated 

from Figure I) 

Included 

Included 

No new data (data from 

O'Reardon et al," 2007) 

No sham (one single sham session 

followed by active treatment) 

DLPFC not stimulated (parietal 

cortex stimulated) 

No new data (data from 

Schutter et al,™ 2009) 

No new data (data from 

O'Reardon et al," 2007) 

Included 

Inadequate data reported (baseline 

scores/group missing) 

Included 

Inadequate data reported (HAMD 
scores missing) 
Included 
Case study 

No sham (sham was the active left- 
slow stimulation of the DLPFC) 
No sham 

No new data (same cases as in 

Zhengetal,"'20IO) 

Included 



Note: A total of N- 1 8/50 studies were included in the final meta-analysis. 
Abbreviations: DLPFC, dorsolateral prefrontal cortex; HAMD, Hamilton 
Depression Rating Scale; rTMS, repetitive transcranial magnetic stimulation: SD, 
standard deviation: 
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Author, year 



George et al,' 1997 
Avery etal,= 1999 
Kimbrell etal,» 1999 all 
Klein etal,' 1999 
Looetal,' 1999 
Padberg etal,= 1999 all 
Berman et al/ 2000 
Eschweiler et al," 2000 all 
George etal,' 2000 all 
Garcla-Toro et al," 2001 
Garcla-Toro et al," 2001 
Manes etal," 2001 
Boutros," 2002 
Padberg etal," 2002 
Fitzgerald et al," 2003 all 
Hoppner et al," 2003 
Loo etal," 2003 
Nahasetal," 2003 
Buchholtzetal," 2004 
Hausmann et al,» 2004 
Holtzhelmeretal," 2004 
Kauffmann et al.^^ 2004 
Koerselman et al,2' 2004 
Mosimann et al,^"* 2004 
Pouletetal," 2004 
Rossini et al," 2005 
RumI etal," 2005 
Su et al,» 2005 all 
Avery etal," 2006 
Fitzgerald et al,'» 2006 
Garcla-Toro et al," 2006 
Januel et al," 2006 
Anderson et al," 2007 
BortolomasI et al," 2007 
Herwig et al,'= 2007 
Loo etal," 2007 
O'Reardon at al," 2007 
Stern et al," 2007 all 
Bretlau et al,™ 2008 
Moggetal," 2008 
George etal," 2010 
Palllere Martlnot et al," 2010 
PallantI et al," 2010 
Triggs et al,'= 2010 all 
Zheng etal," 2010 
Agulrre et al," 2011 
He etal,'' 2011 
LIngeswaran et al," 201 1 
Ray etal," 2011 
Bakim etal," 2012 
Blumbergeretal,*" 2012 
Fitzgerald etal," 2012 
Huang etal," 2012 
Peng etal," 2012 
Chen etal," 2013 
Hernandez-Ribas et al," 2013 
Spamplnato et al," 2013 
Speeretal," 2013 all 



Outcome Statistics for each study 



Std diff in means and 95% CI 



HAMD 
Combined 
HAMD 
Combined 
Combined 
HAMD 
HAMD 
Combined 
HAMD 
Combined 
Combined 
HAMD 
HAMD 
Combined 
Combined 
Combined 
Combined 
HAMD 
HAMD 
Combined 
Combined 
HAMD 
HAMD 
Combined 
Combined 
HAMD 
MADRS 
Combined 
Combined 
Combined 
HAMD 
HAMD 
MADRS 
Combined 
Combined 
Combined 
Combined 
HAMD 
HAMD 
Combined 
Combined 
Combined 
HAMD 
Combined 
Combined 
HAMD 
HAMD 
Combined 
HAMD 
Combined 
HAMD 
Combined 
Combined 
Combined 
Combined 
HAMD 
Combined 
HAMD 



Std diff 
in means 

-1.44 
0.50 
0.44 
0.75 
0.25 
0.58 
1.21 
1.35 
0.70 
0.63 
0.44 
0.30 
0.25 
0.99 
0.44 
0.19 
0.12 
0.09 
0.14 
0.34 
0.22 
0.84 
0.15 
0.12 
0.13 
0.71 
1.23 
1.14 
0.67 
0.53 
0.80 
1.58 
0.78 
1.20 
0.17 
0.34 
0.30 
1.45 
0.75 
0.46 
0.11 
0.40 
1.33 
0.19 
1.85 
0.18 
0.44 
0.09 
2.72 
0.93 
0.22 
0.39 
0.95 
2.45 
0.38 
0.89 
2.04 
0.90 
0.62 



Lower 
limit 

-2.73 
-2.23 
-1.75 
-1.24 
-0.67 
-1.57 
-2.17 
-2.62 
-1.48 
-1.31 
-1.29 
-1.18 
-1.20 
-1.92 
-0.99 
-0.70 
-1.03 
-0.73 
-0.95 
-1.02 
-1.24 
-2.03 
-0.70 
-0.95 
-0.77 
-1.13 
-1.87 
-1.95 
-1.16 
-1.11 
-1.71 
-2.46 
-1.60 
-2.20 
-0.54 
-0.98 
-0.53 
-2.16 
-1.35 
-0.99 
-0.40 
-1.10 
-2.01 
-0.81 
-2.67 
-0.86 
-0.88 
-0.93 
-3.58 
-1.79 
-0.43 
-1.03 
-1.51 
-3.42 
-0.55 
-1.79 
-3.07 
-1.78 
-0.77 



Upper 

limit P-value 

-0.16 0.028 

1.23 0.572 
0.86 0.503 

-0.25 0.003 

1.18 0.592 

0.42 0.258 

-0.26 0.013 

-0.08 0.038 

0.08 0.080 

0.05 0.071 

0.41 0.309 

0.58 0.505 

0.70 0.606 

-0.06 0.038 

0.10 0.110 

1.08 0.676 

0.78 0.787 

0.91 0.828 

1.24 0.796 
0.33 0.318 
0.80 0.674 
0.36 0.171 
0.40 0.564 



0.71 
1.04 



0.774 
0.772 



-0.30 0.001 

-0.60 0.000 

-0.33 0.006 

-0.18 0.007 

0.05 0.075 

0.11 0.086 

-0.71 0.000 

0.04 0.063 

-0.19 0.020 

0.19 0.355 

0.30 0.298 

-0.08 0.009 

-0.74 0.000 

-0.15 0.015 

0.07 0.087 

0.18 0.456 

0.31 0.269 

-0.64 0.000 

0.43 0.548 

-1.03 0.000 

0.50 0.605 

0.01 0.054 

0.75 0.836 

-1.86 0.000 

-0.07 0.035 

-0.87 0.505 

0.24 0.227 

-0.40 0.001 

-1.49 0.000 

1.30 0.427 

0.01 0.052 

-1.00 0.000 

-0.01 0.047 

-0.48 0.000 



Std 

residual P-value 

-1.07 0.28 

0.13 0.90 

0.23 0.82 

-0.26 0.79 

1.42 0.16 

0.07 0.94 

-0.94 0.35 

-0.95 0.34 

-0.13 0.90 

-0.01 1.00 

0.31 0.76 

0.54 0.59 

0.59 0.55 

-0.59 0.55 

0.37 0.71 

1.34 0.18 

0.82 0.41 

1.23 0.22 
1.12 0.26 
0.53 0.60 
0.62 0.54 

-0.29 0.77 

0.96 0.34 

0.86 0.39 

1.24 0.21 
-0.20 0.84 
-1.19 0.23 
-0.90 0.37 
-0.10 0.92 

0.19 0.85 

-0.29 0.77 

-1.60 0.11 

-0.27 0.79 

-0.88 0.38 

1.02 0.31 

0.55 0.59 

0.76 0.44 

-1.54 0.12 

-0.25 0.80 

0.34 0.73 

1.20 0.23 

0.42 0.68 

-1.32 0.19 

0.85 0.40 

-2.12 0.03 

0.84 0.40 

0.40 0.69 

0.91 0.36 

-3.53 0.00 

-0.51 0.61 

1.63 0.10 

0.45 0.66 

-0.67 0.50 

-2.89 0.00 

1.61 0.11 

-0.44 0.66 

-2.14 0.03 

-0.45 0.65 



-3.00 



-1.50 



0.00 



1.50 



3.00 



rTMS effective sham effective 



Figure SI Random-effects meta-analysis of N-58 studies with standardized residuals and their P-values. 

Notes: "AH" indicates that rTMS was administered using different properties into different subgroups of patients in a study and the depression scores for such subgroups 
were combined. "Combined" indicates that more than one depression scale was used in a study and the effect sizes according to the multiple scales were combined). 
According to the P-values, 4/ 1 8 "new" studies were classified as outliers: Zheng et al 20 1 0,^" Ray et al 20 II Peng et al 20 1 2,''' and Spampinato et al 20 1 3.^' These studies 
were excluded from all subsequent analyses. 

Abbreviations: CI, confidence interval; HAMD, Hamilton Depression Rating Scale; MADRS, Montgomery Asberg Depression Rating Scale; rTMS, repetitive transcranial 
magnetic stimulation; Std diff, standardized mean difference d; Std, standardized. 
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Group by 
Outliers 



Outlier 
Study 
Overall 



Statistics for each study 

Std diff Standard Lower Upper 
in means error limit limit P-value 



Std diff in means and 95% CI 



-2.26 0.232 -2.72 

-0.51 0.059 -0.63 

-0.62 0.058 -0.73 



-1.81 0.000 
-0.39 0.000 
-0.51 0.000 



♦ 
♦ 



-1.50 0.00 



1.50 3.00 



rTMS effective sham effective 



Figure S2 Forest plot showing the comparison of the N-4 outliers with N-54 studies. 

Notes: Subgroup analysis using the mixed-effects model revealed that the overall mean weighted effect size d was significantly higher in the N— 4 outlier studies (d— — 2.26) 
compared to the N-54 studies (d--0.5 1 ): Q(df I )-58.3, P<0.00 1 . The variability of the weighted effect sizes was high among the N-4 outlier studies (SEM -0.23) compared 
to the N-54 studies (SEM -0.06). 

Abbreviations: CI, confidence interval; rTMS, repetitive transcranial magnetic stimulation; SEM, standard error of mean; Std diff, standardized mean difference d. 
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Author, year 



Statistics for each study 



Sample size 



Std diff in means and 95% CI 







Std diff 


Lower 


Upper 












in means 


limit 


limit 


P-value 


Sham rTMS 


George etal,M997 


HAMD 


-1.44 


2 


■ 


0 028 


5 


7 


Avety etal,M999 


Combined 


-0.50 


-2.23 


1 .23 


0.572 


2 


4 


Kimbreii et al,^ 1999 aii 


HAMD 


-0.44 


-1.75 


0.86 


0.503 


3 


10 


Klein et al* 1999 


Combined 


-0.75 


-1.24 


-0.25 


0.003 


32 


35 


Loo et al,^ 1999 


Combined 


0.25 


-0.67 


1.18 


0.592 


9 


9 


Padberg et al,^ 1999 ali 


HAMD 


-0.58 


-1.57 


0.42 


0.258 


6 


12 


Berman etal/2000 


HAMD 


-1.21 


-2.17 


-0.26 


0.013 


10 


10 


Eschweiler et al,^2000 all 


Combined 


-1.35 


-2.62 


-0.08 


0.038 


5 


7 


George et al,^2000 aii 


HAMD 


-0.70 


-1.48 


0.08 


0.080 


10 


20 


Garcia-Toro et ai;°2001 


Combined 


-0.63 


-1.31 


0.05 


0.071 


18 


17 


Garcia-Toro et al/^ 2001 


Combined 


-0.44 


-1.29 


0.41 


0.309 


11 


11 


Manes et al,'^2001 


HAMD 


-0.30 


-1.18 


0.58 


0.505 


10 


10 


Boutros,"2002 


HAMD 


-0.25 


-1.20 


0.70 


0.606 


7 


11 


Padberg et al,^''2002 


Combined 


-0.99 


-1.92 


-0.06 


0.038 


10 


10 


Fitzgeraid et al,''^2003 all 


Combined 


-0.44 


-0.99 


0.10 


0.1 10 


20 


40 


Hoppner et al,^^2003 


Combined 


0.19 


-0.70 


1.08 


0.676 


10 


10 


Loo et ai;^2003 


Combined 


-0.12 


-1.03 


0.78 


0.787 


10 


9 


Nahas et ai,^^2003 


HAMD 


0.09 


-0.73 


0.91 


0.828 


12 


11 


Buchhoitz et al,''^ 2004 


HAMD 


0.14 


-0.95 


1.24 


0.796 


7 


6 


Hausmann et al,^°2004 


Combined 


-0.34 


-1.02 


0.33 


0.318 


13 


25 


Holtzheimer et al,^^ 2004 


Combined 


-0.22 


-1.24 


0.80 


0.674 


8 


7 


Kauffmann et al,^^ 2004 


HAMD 


-0.84 


-2.03 


0.36 


0.171 


5 


7 


Koerseiman et 31,^^2004 


HAMD 


-0.15 


-0.70 


0.40 


0.584 


25 


26 


Mosimann et al," 2004 


Combined 


-0.12 


-0.95 


0.71 


0.774 


9 


15 


Pouiet et al,^^2004 


Combined 


0.13 


-0.77 


1.04 


0.772 


9 


10 


Rossini et ai,=^2005 


HAMD 


-0.71 


-1.13 


-0.30 


0.001 


47 


49 


Rumi etal,^^2005 


MADRS 


-1.23 


-1.87 


-0.60 


0.000 


24 


22 


Su et ai,^^2005 ail 


Combined 


-1.14 


-1.95 


-0.33 


0.006 


10 


20 


Avery etal,^^2006 


Combined 


-0.67 


-1.16 


-0.18 


0.007 


33 


35 


Fitzgeraid et al,3°2006 


Combined 


-0.53 


-1.11 


0.05 


0.075 


22 


25 


Garcia-Toro et al,^^ 2006 


HAMD 


-0.80 


-1.71 


0.11 


0.086 


10 


10 


Januel et 31,^^2006 


HAMD 


-1.58 


-2.46 


-0.71 


0.000 


16 


11 


Anderson et ai,3^2007 


MADRS 


-0.78 


-1.60 


0.04 


0.063 


14 


11 


Bortoiomasi et al,^" 2007 


Combined 


-1.20 


-2.20 


-0.19 


0.020 


7 


12 


Herwig etal,^^2007 


Combined 


-0.17 


-0.54 


0.19 


0.355 


59 


57 


Loo et ai,^^2007 


Combined 


-0.34 


-0.98 


0.30 


0.298 


19 


19 


O'Reardon et al,^^ 2007 


Combined 


-0.30 


-0.53 


-0.08 


0.009 


146 


155 


Stem et al,^^2007 all 


HAMD 


-1.45 


-2.16 


-0.74 


0.000 


14 


29 


Bretlau et al,^^2008 


HAMD 


-0.75 


-1.35 


-0.15 


0.015 


23 


22 


Mogg et al,''°2008 


Combined 


-0.46 


-0.99 


0.07 


0.087 


29 


28 


George et al,^^ 2010 


Combined 


-0.11 


-0.40 


0.18 


0.456 


94 


88 


Paillere Martinot et al/^ 201 0 


Combined 


-0.40 


-1.10 


0.31 


0.269 


14 


18 


Pallanti et ai,"2010 


HAMD 


-1.33 


-2.01 


-0.64 


0.000 


20 


20 


Triggs et al,^^2010 all 


Combined 


-0.19 


-0.81 


0.43 


0.548 


14 


34 


Aguirre et al,"*^ 2011 


HAMD 


-0.18 


-0.86 


0.50 


0.605 


15 


19 


Heetal,"2011 


HAMD 


-0.44 


-0.88 


0.01 


0.054 


43 


37 


Lingeswaran et al,^^2011 


Combined 


-0.09 


-0.93 


0.75 


0.836 


14 


9 


Bakim et al,'"'2012 


Combined 


-0.93 


-1.79 


-0.07 


0.035 


12 


11 


Blumberger etal,''^2012 


HAMD 


0.22 


-0.43 


0.87 


0.505 


18 


19 


Fitzgeraid etal."2012 


Combined 


-0.39 


-1.03 


0.24 


0.227 


18 


21 


Huang etal.^2012 


Combined 


-0.95 


-1.51 


-0.40 


0.001 


28 


28 


Chenet al,"2013 


Combined 


0.38 


-0.55 


1.30 


0.427 


10 


10 


Hernandez-Ribas et ai,^2013 


HAMD 


-0.89 


-1.79 


0.01 


0.052 


11 


10 


Speer et ai,«^2013aii 


HAMD 


-0.90 


-1.78 


-0.01 


0.047 


8 


16 






-0.51 


-0.63 


-0.39 


0.000 







-1.50 0.00 1.50 3.00 

rTiUIS effective sham effective 



Figure S3 Random -effects meta-analysis of N=54 studies. 

Notes: "AH" indicates that rTMS was administered using different properties into different subgroups of patients in a study and the depression scores for such subgroups 
were combined. "Combined" indicates that more than one depression scale was used in a study and the effect sizes according to the multiple scales were combined). The 
mean number of patients per group was used in the final calculations if patients dropped out throughout the study between baseline and final sessions. 
Abbreviations: CI, confidence interval; HAMD, Hamilton Depression Rating Scale; MADRS, Montgomery Asberg Depression Rating Scale; rTMS, repetitive transcranial 
magnetic stimulation; Std diff, standardized mean difference d. 
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Author, year Outcome Statistics with study removed Std diff in means {95% 

CI) with study removed 

Lower Upper 
Point limit limit P-value 



George et al,^ 1997 


HAMD 


-0.50 


-0.62 


-0.39 


0.000 


J 


Avery et al,^ 1999 


Combined 


-0.51 


-0.63 


-0.39 


0.000 


S" 


Kimbrell et al,^ 1999 all 


HAMD 


-0.51 


-0.63 


-0.39 


0.000 




Klein et al,"* 1999 


Combined 


-0.50 


-0.62 


-0.39 


0.000 




Looetal,' 1999 


Combined 


-0.52 


-0.64 


-0.40 


0.000 


_ 


Padberg et al,® 1999 all 


HAMD 


-0.51 


-0.63 


-0.39 


0.000 


_ 


Berman etal,'2000 


HAMD 


-0.50 


-0.62 


-0.38 


0.000 




Eschweller et al,®2000 all 


Combined 


-0.50 


-0.62 


-0.39 


0.000 


— 


George et al,^2000 all 


HAMD 


-0.51 


-0.63 


-0.39 


0.000 




Garcla-Toroetal,"2001 


Combined 


-0.51 


-0.62 


-0.39 


0.000 




Garcla-Toroetal,"2001 


Combined 


-0.51 


-0.63 


-0.39 


0.000 




Manes et al,^^2001 


HAMD 


-0.51 


-0.63 


-0.40 


0.000 




Boutros,"2002 


HAMD 


-0.51 


-0.63 


-0.40 


0.000 




Padberg et al,"2002 


Combined 


-0.50 


-0.62 


-0.39 


0.000 




Fitzgerald et al,^^2003 all 


Combined 


-0.51 


-0.63 


-0.39 


0.000 




Hoppner et al,^''2003 


Combined 


-0.52 


-0.64 


-0.40 


0.000 




Looetal, "2003 


Combined 


-0.52 


-0.63 


-0.40 


0.000 


___ 


Nahasetal,"2003 


HAMD 


-0.52 


-0.64 


-0.40 


0.000 




Buchholtz et al," 2004 


HAMD 


-0.52 


-0.63 


-0.40 


0.000 


_ 


Hausmann et al,"2004 


Combined 


-0.51 


-0.63 


-0.40 


0.000 


_ 


Holtzhelmer et al,'' 2004 


Combined 


-0.51 


-0.63 


-0.40 


0.000 


— 


Kauffmannetal,2'2004 


HAMD 


-0.51 


-0.63 


-0.39 


0.000 




Koerselman etal,'=2004 


HAMD 


-0.52 


-0.64 


-0.40 


0.000 




Moslmann et al,"2004 


Combined 


-0.52 


-0.63 


-0.40 


0.000 


— 


Pouletetal,2''2004 


Combined 


-0.52 


-0.64 


-0.40 


0.000 




Rossini etal,»2005 


HAMD 


-0.50 


-0.62 


-0.38 


0.000 




Rumi etal,'^2005 


MADRS 


-0.49 


-0.61 


-0.38 


0.000 


_ 


Su et al,»2005 all 


Combined 


-0.50 


-0.62 


-0.38 


0.000 


_ 


Avery et al,'^2006 


Combined 


-0.51 


-0.63 


-0.39 


0.000 


_ 


Fitzgerald et al,'°2006 


Combined 


-0.51 


-0.63 


-0.39 


0.000 


_ 


Garcia-Toro et al," 2006 


HAMD 


-0.50 


-0.62 


-0.39 


0.000 


_ 


Januel etal,"2006 


HAMD 


-0.49 


-0.61 


-0.38 


0.000 


_ 


Anderson etal,»»2007 


MADRS 


-0.51 


-0.62 


-0.39 


0.000 


_ 


Bortolomasi et al,"2007 


Combined 


-0.50 


-0.62 


-0.39 


0.000 


_ 


Hera/ig etal,"2007 


Combined 


-0.52 


-0.64 


-0.40 


0.000 




Looetal," 2007 


Combined 


-0.51 


-0.63 


-0.40 


0.000 


_ 


O'Reardonetal,"2007 


Combined 


-0.52 


-0.64 


-0.40 


0.000 


_ 


Stem etal,"2007 all 


HAMD 


-0.49 


-0.60 


-0.38 


0.000 


— 


Bretlau et al,®^2008 


HAMD 


-0.51 


-0.62 


-0.39 


0.000 


_ 


Mogg etal,"2008 


Combined 


-0.51 


-0.63 


-0.39 


0.000 


_ 


George etal,^5 2010 


Combined 


-0.53 


-0.64 


-0.41 


0.000 


— 


Paillere Martinot et al,'= 201 0 


Combined 


-0.51 


-0.63 


-0.39 


0.000 


— 


Paiianti et al "2010 


HAMD 


-0.49 


-0.60 


-0.38 


0.000 




Trinnc ot al ^5 9nl n all 
1 iigyb Bi di, lu dii 


Combined 












Aguirre et al,"^ 2011 


HAMD 


-0.52 


-0.64 


-0.40 


0.000 




Meet al,"2011 


HAMD 


-0.51 


-0.63 


-0.39 


0.000 




Lingeswaran et al,^^2011 


Combined 


-0.52 


-0.63 


-0.40 


0.000 




Bakimetal.'"2012 


Combined 


-0.50 


-0.62 


-0.39 


0.000 




Blumberger etal,''^2012 


HAMD 


-0.52 


-0.64 


-0.41 


0.000 




Fitzgerald etal,^2012 


Combined 


-0.51 


-0.63 


-0.39 


0.000 




Huang etal,«2012 


Combined 


-0.50 


-0.61 


-0.38 


0.000 




Chenet al,s*2013 


Combined 


-0.52 


-0.64 


-0.40 


0.000 




Hemandez-Ribas et al,^^ 201 3 


HAMD 


-0.50 


-0.62 


-0.39 


0.000 




Speer et al,^^ 201 3 all 


HAMD 


-0.50 


-0.62 


-0.39 


0.000 








-0.51 


-0.63 


-0.39 


0.000 





-1.00 -0.50 0.00 0.50 1.00 

rTMS effective sham effective 

Figure S4 One study removed analysis (N-54 studies). 

Notes: "AH" indicates that rTMS v^as administered using different properties into different subgroups of patients in a study and the depression scores for such subgroups 
were combined. "Combined" indicates that more than one depression scale v^as used in a study and the effect sizes according to the multiple scales v/ere combined). "Point" 
refers to the overall mean v^eighted d of all studies except for the study listed in each row. 

Abbreviations: CI, confidence interval; HAMD, Hamilton Depression Rating Scale; MADRS, Montgomery Asberg Depression Rating Scale; rTMS, repetitive transcranial 
magnetic stimulation; Std diff, standardized mean difference d. 
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Author, year 



Cumulative Statistics 









Lower 


Upper 












Point 


limit 


limit 


P-value 


Sham 


rTMS 


George et al,' 1997 


HAMD 


-1.44 


-2.73 


-0.16 


0.028 


5 


7 


Avery et al,^ 1999 


Combined 


-1.11 


-2.14 


-0.08 


0.035 


7 


11 


Kimbrell et al,' 1999 all 


HAMD 


-0.85 


-1.66 


-0.04 


0.039 


10 


21 


Klein etal/ 1999 


Combined 


-0.77 


-1.20 


-0.35 


0.000 


42 


56 


Loo etal,= 1999 


Combined 


-0.56 


-1.07 


-0.05 


0.030 


51 


65 


Padberg et al,' 1999 all 


HAMD 


-0.58 


-0.97 


-0.20 


0.003 


57 


77 


Berman etal72000 


HAMD 


-0.67 


-1.04 


-0.29 


0.000 


67 


87 


Eschweiler et al,»2000 all 


Combined 


-0.72 


-1.07 


-0.36 


0.000 


72 


94 


George et al.»2000 all 


HAMD 


-0.71 


-1.01 


-0.41 


0.000 


82 


114 


Garcia-Toro et alj°2001 


Combined 


-0.70 


-0.97 


-0.42 


0.000 


100 


131 


Garcia-Tcro et al," 2001 


Combined 


-0.67 


-0.94 


-0.41 


0.000 


111 


142 


Manes et al,"2001 


HAMD 


-0.64 


-0.90 


-0.39 


0.000 


121 


152 


Boutros,"2002 


HAMD 


-0.62 


-0.86 


-0.38 


0.000 


128 


163 


Padberg et al,"2002 


Combined 


-0.64 


-0.88 


-0.41 


0.000 


138 


173 


Fitzgerald et al,'s2003 all 


Combined 


-0.61 


-0.83 


-0.40 


0.000 


158 


213 


Hoppneretal,"'2003 


Combined 


-0.57 


-0.78 


-0.36 


0.000 


168 


223 


Loo etal," 2003 


Combined 


-0.54 


-0.75 


-0.34 


0.000 


178 


232 


Nahasetal,^'2003 


HAMD 


-0.51 


-0.70 


-0.31 


0.000 


190 


243 


Buchholtz et al,^^ 2004 


HAMD 


-0.49 


-0.68 


-0.29 


0.000 


197 


249 


Hausmann et al,^''2004 


Combined 


-0.48 


-0.66 


-0.29 


0.000 


210 


274 


Hoitzheimer et al,^^ 2004 


Combined 


-0.47 


-0.65 


-0.28 


0.000 


218 


281 


i\auTTiTiann ei ai, ^uu^ 


HAMD 














knorcolman ot al 23 9nn4 


HAM D 














Mosimann et al,^* 2004 


Combined 


-0.43 


-0.60 


-0.26 


0.000 


257 


329 


Poulet etal, "2004 


Combined 


-0.41 


-0.58 


-0.25 


0.000 


266 


339 


Rossini et al,"2005 


HAMD 


-0.45 


-0.61 


-0.30 


0.000 


31 3 


388 


Rumi et al,"2005 


MADRS 


-0.49 


-0.65 


-0.33 


0.000 


337 


410 


Su etal, "2005 all 


Combined 


-0.51 


-0.67 


-0.35 


0.000 


347 


430 


Avery etal,=» 2006 


Combined 


-0.52 


-0.68 


-0.37 


0.000 


380 


465 


Fitzgerald et al,"'2006 


Combined 


-0.52 


-0.67 


-0.38 


0.000 


402 


490 


Garcia-Toro et al,^^ 2006 


HAMD 


-0.53 


-0.67 


-0.39 


0.000 


412 


500 


Januel et al,»2006 


HAMD 


-0.55 


-0.71 


-0.40 


0.000 


428 


511 


Anderson etal," 2007 


MADRS 


-0.56 


-0.71 


-0.41 


0.000 


442 


522 


Bortolomasi etal, "2007 


Combined 


-0.57 


-0.72 


-0.43 


0.000 


449 


534 


Herwig etal," 2007 


Combined 


-0.54 


-0.69 


-0.40 


0.000 


508 


591 


Loo etal," 2007 


Combined 


-0.53 


-0.67 


-0.39 


0.000 


527 


610 


O'Reardon et al,"2007 


Combined 


-0.51 


-0.64 


-0.38 


0.000 


673 


765 


Stem etal," 2007 all 


HAMD 


-0.54 


-0.68 


-0.40 


0.000 


687 


794 


Bretlau et al,"200B 


HAMD 


-0.55 


-0.69 


-0.41 


0.000 


710 


816 


Mogg etal,'°2008 


Combined 


-0.54 


-0.68 


-0.41 


0.000 


739 


844 


George et al," 2010 


Combined 


-0.52 


-0.65 


-0.39 


0.000 


833 


932 


Paillere Martinot et al," 201 0 


Combined 


-0.52 


-0.65 


-0.39 


0.000 


847 


950 


Pallanti etal,"2010 


HAMD 


-0.54 


-0.67 


-0.41 


0.000 


867 


970 


Triggs et al,"2010 all 


Combined 


-0.53 


-0.66 


-0.40 


0.000 


881 


1004 


Aguirre et ai,**^ 201 1 


HAMD 


-0.52 


-0.65 


-0.39 


0.000 


896 


1023 


He et al,"2011 


HAMD 


-0.52 


-0.64 


-0.39 


0.000 


939 


1060 


Lingesvtfaran et al,"201 1 


Combined 














Bakim et al,^ 2012 


Combined 


-0.52 


-0.64 


-0.40 


0.000 


965 


1080 


Blumbergeretal,^'2012 


HAMD 


-0.50 


-0.62 


-0.38 


0.000 


983 


1099 


Fitzgerald etal,'* 2012 


Combined 


-0.50 


-0.62 


-0.38 


0.000 


1001 


1120 


Huang etal," 201 2 


Combined 


-0.51 


-0.63 


-0.39 


0.000 


1029 


1148 


Chen et al,"2013 


Combined 


-0.50 


-0.62 


-0.38 


0.000 


1039 


1158 


Hernandez-Rlbas et al," 201 3 


HAMD 


-0.50 


-0.62 


-0.39 


0.000 


1050 


1168 


Speeretal,"2013all 


HAMD 


-0.51 


-0.63 


-0.39 


0.000 


1058 


1184 






-0.51 


-0.63 


-0.39 


0.000 







Cumulative std diff in means (95% CI) 



-1.50 0.00 1.50 1.00 

rTMS effective sham effective 



Figure S5 Cumulative meta-analysis (N-54 studies). 

Notes: "AH" indicates that rTMS v^as administered using different properties into different subgroups of patients in a study and the depression scores for such subgroups 
were combined. "Combined" indicates that more than one depression scale v/as used in a study and the effect sizes according to the multiple scales were combined). "Point" 
refers to the overall mean weighted d of all studies before and including the study listed in each row. 

Abbreviations: CI, confidence interval; HAMD, Hamilton Depression Rating Scale; MADRS, Montgomery Asberg Depression Rating Scale; rTMS, repetitive transcranial 
magnetic stimulation; Std diff, standardized mean difference d. 
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Regression of total stimuli on std diff in means 
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E 
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0.40 
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-0.56 
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-1.04 
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-2.00 
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-4,680.00 2,376.00 9,432.00 16,488.00 23,544.00 30,600.00 37,656.00 44,712.00 51,768.00 58,824.00 65,880.00 

Total stiuli 
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0) 
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tf) 
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o 



Regression of total stimuli on std diff in means 




-2.00 



-1,680.00 1,776.00 5,232.00 



i,688.00 12,144.00 15,600.00 19,056.00 22,512.00 25,968.00 29,424.00 32,880.00 

Total stimuli 



Figure S6 Univariate meta-regression of the total stimuli/study on the effect size d weighted according to the random-effects model. 

Notes: Total stimuli - stimuli/session x total number of sessions. (A) The top figure shows the results of a significant meta-regression (slope P-0.015) in N-33 studies. 
However, the significance of this regression was due to two largest RCTs by O'Reardon et aP^ and George et al" (depicted as two largest circles on the right-hand side of (A). 
(B) These two RCTs are removed from the analysis in the bottom figure (slope, P-0.208). 
Abbreviations: RCT, randomized controlled trial; std diff, standardized mean difference d. 
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Figure S7 Univariate meta-regression of the stimuli/session on the effect size d weighted according to the random -effects model (two largest RCTs removed). 

Notes: The slope of meta-regression conducted on N-3 I studies remained positive and statistically significant (P-0.018) following the removal of O'Reardon et al" and 

George et al^^ studies. 

Abbreviations: RCT, randomized controlled trial; std diff, standardized mean difference d. 



Table S3 Location (country) where the N=54 studies published 
from 1 997 to August 20 1 3 were conducted 



Rank 


Country 


Number 
of studies 


1 


USA 


16 


2 


Australia 


7 


3 


Spain 


6 


4 


Germany 


4 


5 


People's Republic of China, France, Italy 


3 each 


6 


Austria, Canada, Denmark, UK 


2 each 


7 


Brazil, India, the Netherlands, Switzerland, 


1 each 




Taiwan, Turkey 





Note: N does not add up to 54 because some studies were conducted in more 
than one country. 
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